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Abstract 



Nature-inspired algorithms such as Particle Swarm Optimization and Fire- 
fly Algorithm are among the most powerful algorithms for optimization. In 
this paper, we intend to formulate a new metaheuristic algorithm by combin- 
ing Levy flights with the search strategy via the Firefly Algorithm. Numerical 
studies and results suggest that the proposed Levy-flight firefly algorithm is 
superior to existing metaheuristic algorithms. Finally implications for further 
research and wider applications will be discussed. 

■ Citation detail: X.-S. Yang, "Firefly algorithm, Levy flights and global op- 

timization" , in: Research and Development in Intelligent Systems XXVI (Eds 
M. Bramer, R. Ellis, M. Petridis), Springer London, pp. 209-218 (2010). 

1 Introduction 

o ' 

Nature- inspired metaheuristic algorithms are becoming powerful in solving mod- 
ern global optimization problems [2, 3, 5, 7, 9, 18, 17], especially for the NP- 
hard optimization such as the travelling salesman problem. For example, par- 
ticle swarm optimization (PSO) was developed by Kennedy and Eberhart in 
1995 [8, 9], based on the swarm behaviour such as fish and bird schooling in 
nature. It has now been applied to find solutions for many optimization ap- 
plications. Another example is the Firefly Algorithm developed by the author 
[18] which has demonstrated promising superiority over many other algorithms. 
The search strategies in these multi-agent algorithms are controlled random- 
ization, efficient local search and selection of the best solutions. However, the 
randomization typically uses uniform distribution or Gaussian distribution. 

On the other hand, various studies have shown that flight behaviour of many 
animals and insects has demonstrated the typical characteristics of Levy flights 
[4, 13, 11, 12]. A recent study by Reynolds and Frye shows that fruit flies 
or Drosophila melanogaster, explore their landscape using a series of straight 
flight paths punctuated by a sudden 9O turn, leading to a Levy-flight-style 
intermittent scale free search pattern. Studies on human behaviour such as the 
Ju/'hoansi hunter-gatherer foraging patterns also show the typical feature of 
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Levy flights. Even light can be related to Levy flights [1]. Subsequently, such 
behaviour has been applied to optimization and optimal search, and preliminary 
results show its promising capability [11, 13, 15, 16]. 

This paper aims to formulate a new Levy-flight Firefly Algorithm (LFA) 
and to provide the comparison study of the LFA with PSO and other relevant 
algorithms. We will first outline the firefly algorithms, then formulate the 
Levy-flight FA and finally give the comparison about the performance of these 
algorithms. The LFA optimization seems more promising than particle swarm 
optimization in the sense that LFA converges more quickly and deals with global 
optimization more naturally. In addition, particle swarm optimization is just a 
special class of the LFA as we will demonstrate this in this paper. 

2 Firefly Algorithm 

2.1 Behaviour of Fireflies 

The flashing light of fireflies is an amazing sight in the summer sky in the 
tropical and temperate regions. There are about two thousand firefly species, 
and most fireflies produce short and rhythmic flashes. The pattern of flashes is 
often unique for a particular species. The flashing light is produced by a process 
of bioluminescence, and the true functions of such signaling systems are still 
debating. However, two fundamental functions of such flashes are to attract 
mating partners (communication), and to attract potential prey. In addition, 
flashing may also serve as a protective warning mechanism. The rhythmic flash, 
the rate of flashing and the amount of time form part of the signal system that 
brings both sexes together. Females respond to a male's unique pattern of 
flashing in the same species, while in some species such as photuris, female 
fireflies can mimic the mating flashing pattern of other species so as to lure and 
eat the male fireflies who may mistake the flashes as a potential suitable mate. 

The flashing light can be formulated in such a way that it is associated with 
the objective function to be optimized, which makes it possible to formulate new 
optimization algorithms. In the rest of this paper, we will first outline the basic 
formulation of the Firefly Algorithm (FA) and then discuss the implementation 
as well as analysis in detail. 

2.2 Firefly Algorithm 

Now we can idealize some of the flashing characteristics of fireflies so as to 
develop firefly- inspired algorithms. For simplicity in describing our Firefly Al- 
gorithm (FA), we now use the following three idealized rules: 1) all fireflies 
are unisex so that one firefly will be attracted to other fireflies regardless of 
their sex; 2) Attractiveness is proportional to their brightness, thus for any two 
flashing fireflies, the less brighter one will move towards the brighter one. The 
attractiveness is proportional to the brightness and they both decrease as their 
distance increases. If there is no brighter one than a particular firefly, it will 
move randomly; 3) The brightness of a firefly is affected or determined by the 
landscape of the objective function. For a maximization problem, the bright- 
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ness can simply be proportional to the value of the objective function. Other 
forms of brightness can be defined in a similar way to the fitness function in 
genetic algorithms or the bacterial foraging algorithm (BFA) [6, 10]. 

In the firefly algorithm, there are two important issues: the variation of light 
intensity and formulation of the attractiveness. For simplicity, we can always 
assume that the attractiveness of a firefly is determined by its brightness which 
in turn is associated with the encoded objective function. 

In the simplest case for maximum optimization problems, the brightness / 
of a firefly at a particular location x can be chosen as 7(x) oc /(x). However, 
the attractiveness (3 is relative, it should be seen in the eyes of the beholder or 
judged by the other fireflies. Thus, it will vary with the distance between 
firefly i and firefly j. In addition, light intensity decreases with the distance 
from its source, and light is also absorbed in the media, so we should allow 
the attractiveness to vary with the degree of absorption. In the simplest form, 
the light intensity I(r) varies according to the inverse square law I(r) = 
where I s is the intensity at the source. For a given medium with a fixed light 
absorption coefficient 7, the light intensity / varies with the distance r. That 
is 

/ = I e^ r , (1) 

where Iq is the original light intensity. 

As a firefly's attractiveness is proportional to the light intensity seen by 
adjacent fireflies, we can now define the attractiveness /3 of a firefly by 

/3 = /3 e-^ 2 , (2) 

where /3q is the attractiveness at r = 0. 



3 Levy-Flight Firefly Algorithm 

If we combine the three idealized rules with the characteristics of Levy flights, 
we can formulate a new Levy-flight Firefly Algorithm (LFA) which can be 
summarized as the pseudo code shown in Fig. 1. 

In the implementation, the actual form of attractiveness function j3{r) can 
be any monotonically decreasing functions such as the following generalized 
form 

/3(r) = /3 e-^ m , (m > 1). (3) 

For a fixed 7, the characteristic length becomes T = 7™ 1 / m — > 1 as m — > 00. 
Conversely, for a given length scale T in an optimization problem, the parameter 
7 can be used typical initial value. That is 7 = t^. 

The distance between any two fireflies i and j at Xj and Xj, respectively, is 
the Cartesian distance 



\ J2^k-Xj,k) 2 , (4) 
\ k=l 



where Xj & is the A;th component of the spatial coordinate Xj of ith firefly. For 
other applications such as scheduling, the distance can be time delay or any 
suitable forms. 
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Levy-Flight Firefly Algorithm 

begin 

Objective function /(x), x = (x±, ...,Xd) T 
Generate initial population of fireflies Xj {i = 1, 2, n) 
Light intensity Ii at Xj is determined by /(xj) 
Define light absorption coefficient 7 
while (t <MaxGeneration) 
for i = 1 : n all n fireflies 
for j = 1 : i all n fireflies 
if (I, > h) 

Move firefly i towards j in d-dimension via Levy flights 
end if 

Attractiveness varies with distance r via exp[— jr] 
Evaluate new solutions and update light intensity 
end for j 
end for i 

Rank the fireflies and find the current best 
end while 

Postprocess results and visualization 
end 

Figure 1: Pseudo code of the Levy-Flight Firefly Algorithm (LFA). 



The movement of a firefly i is attracted to another more attractive (brighter) 
firefly j is determined by 

Xj = Xj + /?oe~ 7r ?> (xj - Xj) + a sign[rand - -] © Levy, (5) 

where the second term is due to the attraction while the third term is ran- 
domization via Levy flights with a being the randomization parameter. The 
product © means entrywise multiplications. The sign[rand-|] where rand 6 [0, 1] 
essentially provides a random sign or direction while the random step length is 
drawn from a Levy distribution 

Levy ~ u = t~ x , (1 < A < 3), (6) 

which has an infinite variance with an infinite mean. Here the steps of fire- 
fly motion is essentially a random walk process with a power-law step-length 
distribution with a heavy tail. 

3.1 Choice of Parameters 

For most cases in our implementation, we can take (3q = 1, a £ [0, 1], 7 = 1, and 
A = 1.5. In addition, if the scales vary significantly in different dimensions such 
as — 10 5 to 10 5 in one dimension while, say, —0.001 to 0.01 along the other, it is 
a good idea to replace a by aS^ where the scaling parameters Sk(k = 1, d) 
in the d dimensions should be determined by the actual scales of the problem 
of interest. 
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The parameter 7 now characterizes the variation of the attractiveness, and 
its value is crucially important in determining the speed of the convergence and 
how the FA algorithm behaves. In theory, 7 € [0, 00), but in practice, 7 = 0(1) 
is determined by the characteristic length T of the system to be optimized. 
Thus, in most applications, it typically varies from 0.01 to 100. 

3.2 Asymptotic Cases 

There are two important limiting cases when 7 — > and 7—7-00. For 7 — >■ 0, 
the attractiveness is constant /3 = (5q and r — >■ 00, this is equivalent to say 
that the light intensity does not decrease in an idealized sky. Thus, a flashing 
firefly can be seen anywhere in the domain. Thus, a single (usually global) 
optimum can easily be reached. This corresponds to a special case of particle 
swarm optimization (PSO) discussed earlier. Subsequently, the efficiency of 
this special case is the same as that of PSO. 

On the other hand, the limiting case 7 —> 00 leads to V — > and f3(r) — > 5(r) 
(the Dirac delta function) , which means that the attractiveness is almost zero in 
the sight of other fireflies or the fireflies are short-sighted. This is equivalent to 
the case where the fireflies fly in a very foggy region randomly. No other fireflies 
can be seen, and each firefly roams in a completely random way. Therefore, 
this corresponds to the completely random search method. 

As the Levy-flight firefly algorithm is usually in somewhere between these 
two extremes, it is possible to adjust the parameters 7, A and a so that it can 
outperform both the random search and PSO. In fact, LFA can find the global 
optima as well as all the local optima simultaneously in a very effective manner. 



4 Simulations and Results 
4.1 Validation 

In order to validate the proposed algorithm, we have implemented it in Matlab. 
In our simulations, the values of the parameters are a = 0.2, 7 = 1, A = 1.5, 
and 00 = 1. As an example, we now use the LFA to find the global optimum 
of the Ackley function 



/(x) 



-20exp 
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d I 

- 4\ - exp[- ]T cos(27rx i )] + 20 + e, (7) 



which has a global minimum /* = at (0,0,..., 0). The 2D Ackley function 
is shown in Fig. 2, and this global minimum can be found after about 200 
evaluations for 40 fireflies after 5 iterations as shown in Fig. 3. 

Now let us use the LFA to find the optima of some tougher test functions. 
For example, the author introduced a forest function [19] 



/( x ) = (X) eXP [ -J2 S ' m ( X< i) ' _27F < x i< 27r i 



(8) 



i=l 



5 



2 -2 



Figure 2: Ackley function for two independent variables with a global minimum /* = 
at (0,0). 

which has a global minimum /* = at (0, 0, ...,0). The 2D Yang's forest 
function is shown in Fig. 4. However, an important feature of this test function 
is that it is non-smooth and its derivative is not well defined at the optima 
(0, 0, 0) as shown in Fig. 5. 

4.2 Comparison of LFA with PSO and GA 

Various studies show that PSO algorithms can outperform genetic algorithms 
(GA) [7] and other conventional algorithms for solving many optimization prob- 
lems. This is partially due to that fact that the broadcasting ability of the cur- 
rent best estimates gives better and quicker convergence towards the optimal- 
ity A general framework for evaluating statistical performance of evolutionary 
algorithms has been discussed in detail by Shilane et al. [14]. Various test 
functions for optimization algorithms have been developed over many years, 
and a relatively comprehensive review of these test functions can be found in 

[2]- 

Now we will compare the LFA with PSO, and genetic algorithms for various 
standard test functions. We will use the same population size of n = 40 for 
all algorithms in all our simulations. The PSO used is the standard version 
without any inertia function, while the implemented genetic algorithm has a 
mutation probability of 0.05 and a crossover probability of 0.95 without use 
of elitism. After implementing these algorithms using Matlab, we have carried 
out extensive simulations and each algorithm has been run at least 100 times 
so as to carry out meaningful statistical analysis. The algorithms stop when 
the variations of function values are less than a given tolerance e < 10~ 5 . 
The results are summarized in the following table (see Table 1) where the 
global optima are reached. The numbers are in the format: average number of 
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Figure 3: The initial locations of the 40 fireflies (left) and their locations after 5 
iterations (right). 

evaluations (success rate), so 6922 ±537(98%) means that the average number 
(mean) of function evaluations is 6922 with a standard deviation of 537. The 
success rate of finding the global optima for this algorithm is 98%. 

Table 1: Comparison of algorithm performance 





Michalewicz (d=16) 
Rosenbrock ((2=16) 
De Jong (d=256) 
Schwefel (d=128) 
Ackley (d=128) 
Rastrigin 

Easom 
Griewank 
Yang 

Shubert (18 minima) 



89325 ± 7914(95%) 
55723 ±8901(90%) 
25412 ± 1237(100%) 
227329 ± 7572(95%) 
32720 ± 3327(90%) 
110523 ±5199(77%) 
19239 ± 3307(92%) 
70925 ± 7652(90%) 
37079 ± 8920(88%) 
54077 ± 4997(89%) 



PSO 
6922 ± 537(98%) 
32756 ± 5325(98%) 
17040 ± 1123(100%) 
14522 ± 1275(97%) 
23407 ± 4325(92%) 
79491 ± 3715(90%) 
17273 ± 2929(90%) 
55970 ± 4223(92%) 
19725 ± 3204(98%) 
23992 ± 3755(92%) 



LFA 

2889 ± 719(100%) 
6040 ± 535(100%) 
5657 ± 730(100%) 
7923 ± 524(100%) 
4392 ± 2710(100%) 
12075 ± 3750(100%) 
6082 ± 1690(100%) 
10790 ± 2977(100%) 
5152 ± 2493(100%) 
9925 ± 2504(100%) 



We can see that the LFA is much more efficient in finding the global optima 
with higher success rates. Each function evaluation is virtually instantaneous 
on modern personal computer. For example, the computing time for 10,000 
evaluations on a 3GHz desktop is about 5 seconds. Even with graphics for 
displaying the locations of the particles and fireflies, it usually takes less than a 
few minutes. Furthermore, we have used various values of the population size n 
or the number of fireflies. We found that for most problems n = 15 to 50 would 
be sufficient. For tougher problems, larger n can be used, though excessively 
large n should not be used unless there is no better alternative, as it is more 
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Figure 4: Yang's forest function for two independent variables with a global minimum 
/* = at (0,0). 
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Figure 5: Non-smoothness of Yang's forest function near the global minimum (0,0). 
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computationally extensive. 



5 Conclusions 

In this paper, we have formulated a new Levy-flight firefly algorithm and anal- 
ysed its similarities and differences with particle swarm optimization. We then 
implemented and compared these algorithms. Our simulation results for finding 
the global optima of various test functions suggest that particle swarm often 
outperforms traditional algorithms such as genetic algorithms, while LFA is su- 
perior to both PSO and GA in terms of both efficiency and success rate. This 
implies that LFA is potentially more powerful in solving NP-hard problems 
which will be investigated further in future studies. 

The basic Levy-flight firefly algorithm is very efficient. A further improve- 
ment on the convergence of the algorithm is to carry out sensitivity studies by 
varying various parameters such as (3q, 7, a and more interestingly A. These 
could form important topics for further research. In addition, further studies 
on the application of FLA in combination with other algorithms may form an 
exciting area for further research in optimization. 
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